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In de Rham-Gabadadze-Tolley (dRGT) massive gravity and bi-gravity, a non-minimal matter 
coupling involving both metrics generically re-introduces the Boulware-Deser (BD) ghost. A non- 
minimal matter coupling via a simple, yet specific composite metric has been proposed, which 
eliminates the BD ghost below the strong coupling scale. Working explicitly in the metric formu¬ 
lation and for arbitrary spacetime dimensions, we show that this composite metric is the unique 
consistent non-minimal matter coupling below the strong coupling scale, which emerges out of two 
diagnostics, namely, absence of Ostrogradski ghosts in the decoupling limit and absence of the BD 
ghost from matter quantum loop corrections. 


I. INTRODUCTION 

Developing consistent theories of gravity where the 
graviton has a mass has seen renewed interest especially 
in the last few years [1, 2]. Many studies have been de¬ 
voted to constructing a non-linear extension to the Fierz- 
Pauli mass term [3] while keeping the theory manifestly 
ghost-free at the classical level. The non-linear comple¬ 
tion was accomplished only a few years ago for massive 
gravity [4-7] and then for bigravity [8] (where both met¬ 
rics have their own kinetic term). Consequently, an array 
of phenomenological explorations followed which have fo¬ 
cused on studying cosmological applications of these the¬ 
ories. 

Recently, other questions in the context of these Effec¬ 
tive Field Theory formulations of massive gravity and bi¬ 
gravity have been raised. Since these theories are highly 
non-linear and diffeomorphism invariance is explicitly 
broken, it might not be immediately apparent whether 
the coefficients that determine the interactions are nat¬ 
ural in the t’Hooft sense. In fact, these couplings are 
protected by a quantum analogue of screening sourced 
by the Vainhstein effect [9, 10]. There are reasons to 
believe this is, in fact, a feature of all theories which 
exhibit a Vainshtein mechanism [11] —see, for example, 
Refs. [9, 12-15]. 

The nature of the interactions originated at the quan¬ 
tum level has motivated research into developing consis¬ 
tent couplings to the matter sector [8, 16-20] which would 
not excite the Boulware-Deser (BD) ghost [21-23]. Since 
the phenomenology of these theories is explicitly depen¬ 
dent on the way the massive graviton couples to matter, 
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it is important to detail which couplings are theoretically 
consistent, both classically and quantum mechanically, 
before propagating them towards observational tests. 

While most works have investigated cosmological solu¬ 
tions with a minimal coupling to one of the metrics, it 
is possible that a matter coupling built out of an ad¬ 
mixture of both metrics admits a richer cosmological 
phenomenology. Unfortunately, constructing such a cou¬ 
pling to a matter field generically re-introduces the BD 
ghost [21, 24, 25] already at the classical level. Never¬ 
theless, an exception has been pointed out in the liter¬ 
ature: by combining the two metrics appropriately, one 
can eliminate the BD ghost below the strong coupling 
scale [21, 22[^ —see also Ref. [26[ for a complementary 
derivation in the vielbein language. Ideally, one would 
want a matter coupling where the BD ghost is eliminated 
fully, in which the strong coupling scale is not necessar¬ 
ily a physical cut-off, and thus above the strong coupling 
scale the interactions would be redressed owing to the 
Vainshtein mechanism. It has been argued that this may 
be achieved by switching to the vielbein formulation and 
relaxing the symmetric vielbein condition [23[. However, 
additional considerations have revealed that this is not 
possible for the effective composite metric of Ref. [27[. 

Therefore, given the recent and active investment in 
exploring the ghost-freeness of this matter coupling [21- 
23, 25-28], it is important to understand whether there 
are other non-minimal matter couplings which are also 
ghost-free below or even above the strong coupling scale. 
In this paper, we prove that the non-minimal coupling 
originally proposed in Ref. [21] is the unique composite 
metric that avoids the BD ghost below the strong cou¬ 
pling scale. This has profound consequences on the phe- 


^ In 4 dimensions, the strong coupling scale is the usual A 3 = 
(nP Mp)V3, while in D dimensions it becomes A,d—i = 
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nomenology of these theories [29? -31] since it limits the 
choice of matter couplings in these theories to at most 
two free parameters for each matter sector. 

II. UNIQUENESS OF THE COMPOSITE 
METRIC 

We will be working in D dimensions and consider a 
matter coupling of the form 

^matter — 


A. No Ostrogradski ghosts in the decoupling limit 

We will first impose the condition that the non- 
minimal matter coupling (1) does not give rise to Os¬ 
trogradski instabilities [32] for the scalar mode in the 
decoupling limit. The precise meaning of this limit can 
be found, for example, in Ref. [2]. But for our purposes 
this means that we focus on the scalar Stiickelberg mode, 
taking the limit 

(7) 

fiiu -t with (jf = x°’ - ri‘^^df,TT , (8) 


where are effective, composite metrics for the I-th 
matter sector denoted collectively as . One or both of 
the two metrics g^i, and are assumed to have a stan¬ 
dard kinetic term, corresponding to massive gravity or 
bigravity respectively. We shall consider a strictly local 
composite metric, by which we mean that we only con¬ 
sider point-wise operations, including inverting and 

in constructing (x) out of (x) and fbiAx), 
but not their derivatives or non-local operations. The 
dRGT graviton potential is given by 

D 

G = ^a,c/,(/c), withc/,(/c) = /cf;^/c;]^--/c;:;], ( 2 ) 

where IC^ is defined by 

= (3) 

with the branch choice such that /CJ] —>■ 0 when —>■ 

f^iv can be viewed as the deviation of from g^^i,. 
Without loss of generality, we can choose g^^, and K.^ as 
the elementary building blocks to construct G^fJ. Then 
the most general composite metric is given by 

(4) 

N 

with 

N 

n—O 

where [ ] is the trace of the matrix enclosed, 
p^^jj([/C], [/C^],...) are arbitrary functions of the various 
traces of matrix = JC^IC^ and so on. Note that 

P( 7 V,n) ~ order in K.^. Since = g^p/C]] is 

symmetric in its indices, the indices /i and v in Eq. (5) 
are symmetrized implicitly. 

The strategy of our proof is to impose two consistency 
conditions at different steps to restrict the form of G^fJ 
to that of Ref. [21]: 

= 3pp( 4)^^ + 2a(p)%)/C^ + (/C^)^), (6) 

where «(/) and /?(/) are constant. When /3(/) = 0 or 
/3(/) = — 0 ( 7 ), a minimal matter coupling is reproduced. 


where tt is the scalar Stiickelberg mode. In this limit, we 
have 


/c]] = n;; = 


(9) 


whilst the tensor and vector modes are suppressed. We 
will use to lower the indices in this subsection. Con¬ 
sequently, we will use and If]] interchangeably in this 
subsection. G^^} are now simply functions of (and 
VfJ-v)- 

To avoid the Ostrogradski ghost [32], we shall require 
the contribution to the tt equation of motion coming from 
C-matter to not Contain higher order derivatives, either in 
TT or in the matter fields. Since the contributions from 
different Cm{G\iv contain different matter 

fields, these different contributions do not cancel each 
other in the tt equation of motion. This implies that we 
can focus on one matter sector, and, omit the index / 
here and afterwards. Thus, the contributions to the tt 
equation of motion arising from CmiG^^, ip, dip), 


F — c) c) 

C-TT - G/pU(y 


^Q rjifll' dG^if 


5n 

dG 


pa 


= dpd^iV^ + 2a(p(y^ rnd.) 


pa 


+ rn dpd, 


dG, 


dHo 


dG^, 

dTlp, 

( 10 ) 


should not contain higher order derivatives, where G is 
the determinant of G^i, and the energy momentum tensor 
from the I-th matter sector is given by 


-2 d£^{Gp,,ip,dip) 
yZG dGp, • 

It is usually assumed that the matter sectors are diffeo- 
morphism invariant separately, so we have the energy 
momentum conservation for each sector 


= 0 , ( 12 ) 

where is the Christoffel coefficients associated with 
the metric Gpi,. Now, since contains terms with 
first (or higher) derivatives of the matter fields, the first 
term of Eq. (10) contains terms with third (or higher) 
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derivatives of the matter fields, which cannot be canceled 
by other terms. Therefore, as a necessary condition, we 
impose the first term of Eq. (10) to vanish identically: 


dpdaiV^ TPn^ = 0- 

(13) 

For later convenience, we define 


r^^%a = dpdaiV^ TP^ , 

(14) 

= dp{V^ TP’') , 

(15) 

TP’' = TP’' , and 

(16) 

p, pa dGpn 

~ BHp, • 

(17) 

Since TP" p and TP’' pa are first and second derivatives of 


Tfiv respectively, their numerical values at a specific, ar¬ 
bitrarily chosen point in spacetime would be independent 
from that of in the absence of Eq. (12). That is, by 
choosing the matter configuration appropriately in the 
neighborhood of a specific point, the numerical values of 
pa,T^''p and can be assigned independently at 

that point, subject to the following constraints 

- T^'^dpr^f, (isa) 
= -T^pr.p^^. (18b) 


In other words, Eq. (21) has to be an identity. We will 
make use of this identity to constrain the form of G^v 
Now, since this is an identity, different orders of 11, 
and thus different orders of /C, should cancel separately, 
so it is sufficient to consider the iV-th order terms of the 
general ansatz 

N 

= gppY.P^M, (23) 

n—0 

where pAr_„([/C], [/C^],...) are of order 

As an identity, Eq. (21) should be solved by any 
configuration of tt. To further simplify our discus¬ 
sion, it is sufficient to choose a diagonal configura¬ 
tion for HP. (An alternative point of view, which also 
works for our purposes, is that can always be di¬ 
agonalized via an appropriate coordinate transformation 
around a given point.) Suppose {Ao,Ai,--- ,Ad_i} are 
the diagonal components of nj). Then, the / 3 / 3 ““(a ^ 
/3, no summation for aa and /3,5) component of = 

0 gives 

~ B 

r>N Ota. r>N aa ^aa ^ /^N _ n /0/iA 

ypp =yp /3 =v 


To simplify our discussion, we choose a matter configu¬ 
ration where and vanish at a spacetime point. 

This can always be achieved as follows: Suppose there 
are two matter configurations where T'^'^p\p = p\p 
and T'P''\p = TP'^\p, where p is a spacetime point; That 
is, AT^'^pIp = 0 and AT^'^jp = 0; Then, one takes 
Af^ = —2dpda[^TP''dGp,y/dHp„] as our starting 
After this, at that point the constraint system reduces to 
a condition that is much easier to handle: 

= 0. (19) 

To extract the conditions on Gpi, encoded in Eq. (13), 
we need to project out the traces of pa- That is, we 
need a projector V such that pv vanishes for un¬ 

constrained TP'^pa- This reduces the system of Eqs. (13) 
and (19) to a single equation: 

Gp.^'^(VTr%a = 0- (20) 

One can shift the projector V to act on Gpu^^ instead 
and leave TP''pa to be a generic tensor. Then, getting 
rid of the generic tensor TP"pa, Eq- ( 20 ) reduces to the 
requirement: 


Since 77 ““ = ±1, we have BG^p/dXa = 0. Thus, G^^ 
must be independent of Xa{a ^ (3). Since G^p is N-th 
order in 11 , it must be of the form 

G^p = GpX^p, (25) 


where Cp is a constant. It follows from Lorentz invariance 
that the only possible form of Cp should be Crjpp, with 
G being a constant. Since we have (n'^)O = X^Sp (no 
summation for p) for the diagonal configuration chosen, 
we must have 


G;r. = C77pp(n^)e = c(n^)p.. (26) 

Again, by Lorentz invariance, these relations must be also 
satisfied by the non-diagonal components of HP. 

Therefore, there is only one term, PNgpp{K.^)v, at N- 
th order that survives the consistency check in the de¬ 
coupling limit, and we end up with 

G^^ = E Ep^(^'^)- (27) 

N N 


Gp.P’^ = {VG)pu^^ = 0 . ( 21 ) 


In D dimensions, such a projector is explicitly given by^ 


7 




■ po- _ x'^Pn ^ <7)1 

L> + 2 


where pn now are constant. In summary, we have re¬ 
duced our ansatz in Eq. (5) to the one in Eq. (27) by 
requiring that the non-minimal matter coupling does not 
give rise to Ostragradski ghosts in the decouling limit. 


1. Example: Lowest orders 


^ See Appendix A for details on the derivation. 


Before moving on to the next step of the proof, it is 
instructive to give a concrete example to illustrate how 
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these abstract arguments work in essence. Consider the 
most general composite metric, up to 2nd-order in /C: 

G “t“ + a2K.^y + bi[K.'^]g fll^ 

+ b2[}C]^g^i, + &3[/C]/C^i/ + biJC^pg^'^lCa-v, (28) 


where ai and bi are constants and /C^ is to be evaluated 
in the decoupling limit as n)). The identity in Eq. (21) 
can be straightforwardly calculated 




2a 


- "2 G 


D + l '''> D 

+ 2biUP'^gp^ + 2b2[U]gP'^gpi, + 6377 ^'^n^i, 
46i-462(T> + 2) + 263 


+ 

= 0 . 


{D + l){D + 2) 


4n]<5f 


(m 


(29) 


When D > 2, all the terms in Eq. (29) cannot cancel 
each other, so Eq. (29) enforces ai = 61 = 62 = t >3 = 0. 
Thus, up to 2nd order, the consistency requirement in 
the decoupling limit implies 


Gpi, — gpu + a2JCpi, + b4lCppg^'^ICav ■ (30) 

When D = 2, we can get the same result, but one 
needs to take into account the Cayley-Hamilton theo¬ 
rem when checking cancellations between the terms in 
Eq. (29). (The Cayley-Hamilton theorem states that: 
suppose that p{X) = 0 is the characteristic polynomial of 
matrix A, then substituting A for A in the polynomial 
gives rise to an identity, p{A) = 0. To make use of this 
identity in Eq. (29), one can differentiate the identity 
with respect to A: dp{A)/dA = 0.) When N = D, the 
Cayley-Hamilton identity is used directly; when N > D, 
one multiplies p(A) = 0 with powers of A to get relevant 
identities. The diagonalization of H^ in the last subsec¬ 
tion, on the other hand, is a convenient way to avoid the 
complications due to the Cayley-Hamilton identities for 
N>D. 


B. No BD ghost from matter loop corrections 


to 


Gp, = gpp{a^SP + 2aPKP + . (33) 

First, notice that po should be non-zero (positive defi¬ 
nite if the signature of the metric is taken into account), 
otherwise the effective metric Gp^, becomes singular in 
the limit 0. Therefore, we can re-write 

Gp.=Pegpp{5P + PP{K)) (34) 

= Po 9pp p'n ^ (35) 

where p(y = Pq^pn- Note that p^ are constant here. 
Then the determinant of ansatz (34) can be re-cast as 

= pfv^ detiY^ + PilC)) (36) 

= P(f y^det(l-f Q(/C)) (37) 

= pf (j^UYQilC))^, (38) 

where 

gO(/c) = ^ (39) 

N=1 

and qn can be expressed in terms of p'j^ by Taylor ex¬ 
panding yA + PjjC) and comparing to the coefficients 
of {JC^)P. For the requirement (31) to go through, the 
following equation 

D D 

Y,Us{Q{JC)) = J2<^sUs{JC) (40) 

s—0 s—0 

should be satisfied for some constant as- We will check 
what this requirement implies order by order in 1C. 

The 0-th order equation can be satisfied by setting 
oq = 1. At order 1, we have qiUi{lC) = ai[/i(/C), which 
gives qi = ai- At order 2, we have 


Given an effective composite metric Gpi, in dRGT mas¬ 
sive gravity or bigravity, it is natural to include the cos¬ 
mological term V—GA, A being constant, in the La- 
grangian. If it is not there in the bare Lagrangian, it has 
been shown that matter loop corrections will generically 
generate a cosmological term for the effective metric [ 10 ], 
much like that in general relativity. So, to avoid matter 
quantum corrections to re-introduce the BD ghost, we 
require 

D 

s=0 

where Og are constants and UsiJC) are defined in Eq. (2). 
We will show that this requirement is sufficient to reduce 

=Y.^pu= 9pp E (32) 

N N 


qlU2{K.) + q2Ui{lC^) = a2t/2(/C), (41) 

which leads to 

ql = 02, 92 = 0 . (42) 

At order 3, making use of the fact that g 2 = 0, we have 
9?C/3(/C) + 93t/i(/C3) = a3t/3(/C), (43) 

which leads to 

9i =03, 93 = 0 . (44) 

This can be extended to arbitrary orders, so that, at 
order s, we simply have 

1 < s < D : 9 ?[/g(/C) + qsUi{KY = asU^K.), (45) 
s > D : q^UiilCY = 0. (46) 
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That is, by solving Eq. (40) order by order, we can con¬ 
clude that 

qi = ai, Qs =0 (s > 1) , (47) 

which leads to 

Gf,u = Pogtip [(1 + qiK-f] I ■ (48) 

Therefore, by requiring the BD ghost does not re-emerge 
under matter loop corrections, we have reduced the 
ansatz in Eq. (32) to Eq. (33), as advertised. This is pre¬ 
cisely the effective composite metric initially proposed in 
Ref. [21], and it emerged here naturally from requiring 
absence of Ostrogradski ghosts in the decoupling limit 
and absence of the BD ghost from matter loop correc¬ 
tions. 


III. CONCLUSION 

In this letter we have explored a generic class of com¬ 
posite couplings to matter in dRGT massive gravity and 
bigravity, which involve a generic admixture of the met¬ 
rics Qfii, and /^[^. We have imposed two diagnostic tests 
to ensure the ghost-freeness of the theory at least be¬ 
low the strong coupling scale. First, we have required 
such non-minimal coupling not to give rise to Ostrogradki 
ghosts for the scalar Stiickelberg mode in the decoupling 
limit. This has allowed us to discard a big subset of all 
such couplings. Furthermore, we have imposed that mat¬ 
ter loop corrections do not re-introduce the BD ghost, 
which has allowed us to single out one composite met¬ 
ric with two free parameters as the unique non-minimal 
coupling to matter. This is precisely the composite met¬ 
ric proposed in Ref. [21], and given in the metric lan¬ 
guage in Eq. (33). Consequently, cosmological solutions 
in these theories only depend on a finite choice of healthy 
couplings to matter at energy scales comparable to the 
strong coupling scale. 

We note that our proof does not assume any specific 
form of the matter fields in the non-minimal matter cou¬ 
pling - One can view tp as a, vector field encompassing 
all possible fields. However, we do assume that mat¬ 
ter fields are not coupled to derivatives of the metrics 
and in each matter sector all matter fields couple to one 
universal composite metric. More generically, however, 
one may consider a case where derivatives of the metrics 
also enter the matter couplings and matter fields cou¬ 
ple to the two metrics in a convoluted non-trivial way. 
Whether ghost free non-minimal couplings generically ex¬ 
ist in this case is beyond the scope of this letter, but such 
an exotic theory should exist and be free of the BD ghost 
below the strong coupling scale, based on the results of 
Ref. [23, 27] The non-minimal matter coupling (33) has 
a very simple representation in the vielbein formulation. 
Taking this vielbein non-minimal matter coupling as the 


starting point, one can reproduce the non-minimal cou¬ 
pling in the metric formulation (33), if one imposes the 
symmetric vielbein condition [27]. On the other hand, 
if one imposes a modified vielbein condition, as in [23], 
one would end up with a convoluted metric theory that is 
physically different from the theory with (33). However, 
despite being complicated and exotic, the theory with the 
modified vielbein condition has the same decoupling limit 
as the theory with the symmetric vielbein condition [27]. 

Note added: The uniqueness of the composite metric in 
the vielbein formulation has been argued in Ref. [27], 
which appeared when our paper was being finalized. Our 
proof in the metric formulation is complementary to the 
comments in Ref. [27]. 

ACKNOWLEDGMENTS 

We would like to thank Claudia de Rham, Kurt Hin- 
terbichler, Rachel Rosen and Andrew J. Tolley for useful 
discussions. RHR acknowledges the hospitality of the 
Perimeter Institute of Theoretical Physics during stages 
of this work. RHR’s research was supported by a DOE 
grant de-sc0009946, the STFC grant ST/J001546/1 and 
in part by Perimeter Institute for Theoretical Physics. 
Research at Perimeter Institute is supported by the Gov¬ 
ernment of Canada through Industry Canada and by the 
Province of Ontario through the Ministry of Economic 
Development & Innovation. SYZ acknowledges support 
from DOE grant DE-SC0010600, and would like to thank 
Kavli Institute for Theoretical Physics China at the Chi¬ 
nese Academy of Sciences for hospitality during part of 
this work. QGH, YHX and KCZ are supported by the 
grants from NSFC (grant NO. 11322545 and 11335012). 

Appendix A: Projector V 

Consider a generic (2,2) tensor A]]] where the up and 
down two indices are symmetric respectively. We want 
to derive a projector that projects out the traces in A]][. 
That is, for generic A]]], we need {VA)p'[ = 0. Consider¬ 
ing the index structure of A]]], the only trace terms are 

and Therefore we have 

{rA)[Z = + ciAl\l5:] + bAll5[l6':l , (Al) 

where a and b are constants. To determine a and 6, we 
impose the condition that the right hand side of Eq. (Al) 
vanish identically when a and iz are contracted. This 
leads to, in D dimensions, 

^_4_ 2 

“ D + 2' {D + l){D + 2) ■ 
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